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Abstract 

We have investigated the model dependence of meson resonance properties extracted from the 
Dalitz-plot analysis of the three-pions photoproduction reactions on the nucleon. Within a unitary 
model developed in Phys. Rev. D 84, 114019 (2011), we generate Dalitz-plot distributions as 
data to perform an isobar model fit that is similar to most of the previous analyses of three-pion 
production reactions. It is found that the resonance positions from the two models agree well when 
both fit the data accurately, except for the resonance poles near branch points. The residues of 
the resonant amplitudes extracted from the two models and by the usual Breit-Wigner procedure 
agree well only for the isolated resonances with narrow widths. For overlapping resonances, most 
of the extracted residues could be drastically different. Our results suggest that even with high 
precision data, the resonance extraction should be based on models within which the amplitude 
parametrization is constrained by three-particle unitarity condition. 

PACS numbers: 13. 25.-k,14.40.Rt, 11.80. Jy 
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FIG. 1. (color online) M*-decay amplitude. 



I. INTRODUCTION 

Meson properties are important information for understanding the confinement mecha- 
nism of QCD. Thus the investigation of meson spectroscopy has long been an important 
subject in hadron physics. In recent years, more emphasis has been placed on the study of 
mesons with quantum numbers beyond the classification of conventional constituent quark 
model. Such mesons, called exotic mesons, are expected to have explicit gluonic and/or four- 
quark components in their structure Therefore the search for exotic mesons has been 
an important goal in the experiments on irN — > M*N — > mrirN at BNL 0] and CERN 0], 
and 7iV -+ M*N ->■ ttttttN at JLab H, where the intermediate excited mesons M * could be 
exotic. However, the existence of exotic mesons, such as 7Ti(1600) (J pc = 1~ + ) has not been 
conclusive so far. The forthcoming experiments to be performed at JLab after the 12 GeV 
upgrade [5| are aimed at providing high precision data for making progress in this direction. 
In addition to searching for exotic mesons, the new data can also be useful for investigating 
some mesons that could have exotic structure and can be revealed in their characteristic 
decay patterns, as discussed in Ref. [o| with the 3 Pq model. 

We are here interested in the excited mesons that decay into three light mesons 
ttttK, etc.). Since these excited mesons are unstable and couple with multi- mesons contin- 
uum to form resonances, the meson spectroscopy can be determined only by analyzing the 
resonances extracted from the meson production reaction data. Conventionally, these data 
were analyzed by using the so-called isobar model (IM) in which two of the three mesons 
form a light flavor excited meson R (f , p, K*, etc.) and the third meson is treated as 
a spectator of the propagation and the decay of R into two light mesons. This approach 
obviously violates the three-body unitarity and neglects the coupled- channels effects since 
the outgoing R can have multiple scattering with the third meson, as illustrated in Fig. [TJ 

Recently we developed a unitary coupled-channels model j6|. In this work, we use our 
model to analyze the 7p — > M*n — > 7T + ii + ii^n reaction, and study the importance of three- 
body unitarity and coupled-channels effects in extracting the excited meson properties. We 
consider the reaction where the intermediate M* can be several and can overlap. We will 
show that while the IM can fit the same Dalitz plot data generated from our unitary model 
(UM), the extracted resonance parameters are rather different. Our finding indicates the 
limitation of the IM in establishing the meson spectroscopy. 

This paper is organized as follows. In Sec. [TT1 we present formulae based on our model for 
identifying the meson resonances in the amplitudes of the 'yp — > M*n — > 7r + 7r + 7r~ reaction, 
and give the expression for calculating the corresponding Dalitz plots of the cross sections. 
Our procedure and numerical results are presented in Sec. IHfl followed by a summary in 
Sec. ED 
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FIG. 2. (color online) A graphical representation of the 7iV — > abcN' reaction. The momentum 
for each particle is shown to explain the notation used in this work. The white square indicates 
the dressed vertex more explicitly shown in Figs. [3]and[31 



II. FORMULATION 

In this paper, we denote pions by a, b, c and the light excited mesons (such as /o, p ) by 
R. As illustrated in the left-hand-side of Fig. |2l we assume that the 'yN — > abcN' reaction 
proceeds via the photo-production of resonant M* states which decay into cR states. Under 
the three-particle unitarity condition, the propagating cR states experience the multiple- 
scattering due to the Z-diagram mechanism, as illustrated in Fig. [TJ The final three pion 
states are then generated via the R — > tttt decay mechanism. In the vector meson dominance 
(VDM) model, one of the production mechanisms can be calculated from a irp — > M* and 
the well known ttNN vertex, as illustrated in the right-hand-side of Fig. [21 For simplicity, 
we will calculate the cross section of 'yN — > abcN' reaction using this VDM mechanism, and 
the up —> M* vertex that reproduce partial width predicted by the 3 Pq model of Ref. 0]. 
Accordingly all M* — > ttR couplings included in our calculations are also fixed in the same 
manner. Obviously, this is a simplification, but is sufficient for our present purposes in 
investigating the importance of three-particle unitarity in extracting M* from the three- 
pion distribution data. For analyzing the data from CLAS, we need to develop models of 
other mechanisms. This non-trivial task is beyond the scope of this investigation. 



A. Cross section 

With the momentum variables specified in Fig. [2(b), the cross section of 7A^ — > ttttttN 
can be written as 

v rc \ 2E~ 2En 4 ^ 1 



ttttttN' I 



(o \4x(4)( _ ^ d3pN ' d2pa d3pb d3pc m 

x { zn) ^ p f) (27r)32 ^ (27r)32 £ a (2n) s 2E b (2tt) 3 2£ c ' 1 } 

where v re \ is the relative velocity of the initial particles. We used and will use a notation 
E x = ml + pi to denote the free energy for a particle x with the mass m x and the 
momentum p x . We denote the Bose factor for the final three pions by B, and B — 1/2 
for 7r + 7r + 7r~ final states. We have taken the average of the initial nucleon and photon 
polarizations (| ^2 SN x ) an d the summation over the final nucleon spin (%2 8 ,)■ 

We perform the calculations in the center of mass (CM) frame of the total system. The 
orientation of the three final pions are specified by the Euler angles fi Euler = (a, f3, 7) in the 
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FIG. 3. (color online) A graphical representation of the dressed M* — > cR vertex 
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FIG. 4. A graphical representation of cR scattering equation. 

three-meson CM frame (see Eq. (10) and Fig. 1 of Ref. [§[ for the definition of the Euler 
angles). With some manipulation of the four-body phase-space factor Eq. flU) leads to 



da 



Bm 2 N 



dt dW dfW dm\ h dm 2 bc 4 ^{ Pl -p N ) 2 (4tt) 7 E tot p™W 



'yN—HTinrN' | ) 



(2) 



! JVi s jv')-^7 



where <if2 Euler = da <icos/3 dj, E tot is the total energy, p™ is the initial nucleon momentum, 

t = (pn-Pn') 2 , W = v 7 (Pa+Pb + Pc) 2 , and m|- = (Pi+Pj) 2 . For a given set of (t, PF, fi Eu ier), 
we then can calculate the Dalitz-plot distributions of the three out-going pions as functions 
of two-particle invariant masses m ab and m bc by using Eq. fl5]) • 



B. Matrix element 

With the normalization (k\k ') = 5(k — k') for the plane-wave state, the invariant ampli- 
tude A^ 7 Ar^ 7r7rvr Ar' in Eq. ([2]) is related to the scattering amplitude T abcN >^ N by 

M lN ^abcN' — ~TZ — 7^7 

x v /(27r) 3 2E a (p a ) v /(27r) 3 2£; b (p fe ) v /(27r) 3 2E c (p c ) v /(27r) 3 2^(^') T abcN ,^ N (3) 

With the VDM process illustrated in the right-hand-side of Fig. |2} the scattering amplitude 
in Eq. ([3]) can be written within our model as 

cyclic 

T abcN ,,, N = EEE E F{a>vy,M ta {W) [G M *{W)] iaja T M ; aN ,„ N {W) , (4) 

a i a j a S^j* (a'b'c') 

where T M * N > )jN is the M* photo-production amplitude, F( ab y iM * (W) is the M* — > abc decay 
amplitude, and Gm*(W) is the propagator of the M* resonant state. The index a labels a set 
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of quantum numbers (spin, isospin, parity) of M*, and i a ,j a runs over all M*s belonging to 
the set of quantum numbers a. The summation over the M* spin orientation is denoted by 
^2 S z t ■ The symbol ^(a' C &v) means taking summation over the cyclic permutation, (a'b'd) = 
(abc), (cab), (bca). As illustrated in Fig. [31 the M* decay amplitude in Eq. OH) consists of 
two terms 

F{ab)c,M* (E) = Fffi c>M . (E) + F™ c Mt (E) . (5) 
The direct decay amplitude [Fig. 13(a)] is defined by 

F($c,m*( E ) =^^<ab\fabjtGcM«(E)r &RftM .\M*), (6) 

R c'R' 

where f a b,R is the R — > ab vertex interaction, and T^r^m* \M*) describes the M* — > c'R' 
decay. The summation in Eq. ([6]) runs over the particle species and its momentum, spin, 
and isospin components. The second term of Eq. dSJ) [Fig. Mb)] includes the final state 
interaction (FSI), as required by the three-particle unitarity condition. It has the following 
expression: 

r (ab)c,M* 

( S ) = E E E (ab\f ah>R G cR , dR ,{E)T' c IR , c m RIII {E)G C >» R m iC n R n (E)Y ' c » R „ >M * \ M*) , 
R c'R' c"'R"',c"R" 

(7) 

where T' C , R , cR (E) is the cR — > c'R' scattering amplitudes. As illustrated in Fig. HI T' C , R , cR (E) 
is defined by the following coupled- channels scattering equation, 

T'c!R',cr{ e ) — Z c'R',cr(E) + Z dR , :C m R ,,,(E)G c ,,, R » l ^, R ,i(E)T' cllR ,, cR (E), (8) 

c"'R"',c"R" 

where Z c r R > tCR is the one-particle-exchange Z-diagram interaction that is also determined by 
the R — > ab vertex interaction f a b,R 

Z c 'R',cr(E) = y~]fR> >a J>— —fc'c»,R , (9) 

* — ' a, — hi c — hi c i — hi c ii + te 

c" 

where c" is the exchanged meson. 

The Green function in Eqs. (JHJ)- (Ej) is defined by 

[G-\E)] clR , cR = 8 d , c [{E -E c - E R )6 R ,, R - V R ,, R {E - E c )} , (10) 

where the self-energy of light excited-meson R is determined by the R ab vertex interac- 
tion 

sw«o= E — P Bab P ■• M r\ r )- ( n ) 

— _ w — h a — hi), + it 

ab=inr,KK 

Here B a b is a factor associated with the Bose symmetry of mesons: B a b = 1/2 if a and b are 
the identical particles or otherwise B a b = 1. 
The M* propagator in Eq. fl4]) is defined by 

[G-^{E)].. = {E- M° M: )5 tJ - [E^(J5)] y , (12) 
5 



where M™* is a bare mass, and the self energy is given by 

i 

[E M *(E)]ij = ^2 (M*\T M ^ cR G cRfi i R i(E)\T c i R i )M *). 



(13) 



cR,c'R' 



The self energy [T,M*{E)]ij is non- vanishing only when M* and M* have the same spin, 
isospin and parity In Eq. ( fl3l) . \T c r s m*) is the dressed M* — y cR vertex function defined by 



cR,M* 



c'R' 



cR.c'R' 



^2 T c R,c"R"Gc" 



R",c'R 



(E) 



c"R" 



rv 



R',M* 



\M*). 



(14) 



Note that the M* —y abc amplitude F^^m* defined by Eqs. (j3])-([7j) is identical to the matrix 

element J2r c >r> { ab \fab,RG C R, c 'R'\^c'R',M*)- 

To proceed, we need to define R <-y ab and M* — y cR vertex functions. The R — y ab 
vertex functions f a b,R have been determined within our model j(| by fitting the irir scattering 
amplitudes up to 2 GeV with R = /o, p, f2- We include nir and KK channels in our model of 
7T7T scattering. While the tttttt three-body states appear in Z-diagram interactions in Eq. fl9]), 
Ti KK states appear only in the Green function defined by Eq. (fTUl) . 

We thus will only consider cR = nfo, np, tt f'2 channels. For the bare M* <-y cR interaction 
F c r,m*(Pc), we use a partial wave expansion, 

TcrMPc) = ^(t c t^|T M ,T^)(// 2 , s R s z R \S M ^S^)Y* lz (-p c )r icR)ltM ,(p c ) , (15) 



where the first (second) parenthesis is isospin (angular momentum) Clebsch-Gordan coef- 
ficient; s x (t x ) is the spin (isospin) of a particle x and s z x {t z x ) is its ^-component, and I is 
the orbital angular momentum of the relative cR motion. The vertex function Tr cR ) h M*(p) 
is given with the following parametrization for each partial-wave: 



(cR),,M 



(P) 



1 



(2n)^ C(cR) " M ^ AE R (p)E n (p) 



m R 



A 2 



2+(Z/2) 



A 2 



V 



ir\ ■ (16) 



The parameters G( c r\ u m* an d ^{ c r) u m* are chosen to reproduce the partial decay widths of 



M* — y tiR predicted by the 3 Po model of Ref . [9| , and will be explained in the next section. 

To calculate the 7iV — y nnixN' amplitude defined by Eq. (jl]), we also need to calculate 
the VDM photo-production amplitude Tm'N'^n illustrated in the right-hand-side of Fig. |2j 
We use the following interaction Lagrangian to describe the emission of the pion from the 
nucleon: 



NN 



(17) 



where ft NN (q 2 = 0)/4tt = 0.08 and W(g 2 )/W(0) = l/(l-q 2 /Al NN ) 2 .We use A nNN = 
700 MeV which is close to most of the values from the irN scattering model [13(. The photon- 
p contact interaction within the VDM model is defined by the following Lagrangian: 



C 



em 



pi 



9 P 



p p\A» 



(18) 



with g 2 J(4n) = 2.2. 
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With the Lagrangians Eqs. (fT7|) and (|T8|) . the '-/N — > M*N' production amplitude in ^yN 
CM frame is of the following form 



Tm<n'„n(pm»Pn'',P-,,Pn) = \ OF T-i 7^ \ 7-1 1-=, \ y^\J^E Pi (q)E v (q)T M *^ P M ~ 

1 ifnNN(q 2 ) 



q 2 — ml m Tl 



fa) Up N ,i^up N , (19) 



where q = p^ — Pn>, up is the nucleon spinor. The summation is taken over the first and 
second bare p states. The isospin matrix element between the nucleon states is denoted by 
(t x ) = 1 (V2) for N = N' = p (N = p, N' = n). The M* p.n vertex function T M * t7rp .(q) 
has been defined by Eqs. f lT5|) and f|T6|) . 



C. Isobar-model 



We can obtain a model similar to the commonly used IM from the above formula by 
neglecting the nR final state interactions, FFJl)cM* m ^1- ©' which are due to the Z- 
diagram mechanism, as illustrated in Fig. [3](b) and Fig. HI Furthermore, the M* propagator 
Gm*{W) in Eq. 03]) are replaced by a Breit-Wigner (BW) form. Explicitly we consider the 
following expression of the IM, 

cyclic 

TZn',,n = EEEE F ™y,M; a (w) [G B J!(W)} jaja T M;aN ,, 7N (w) , (20) 

a j a S* M * (a'b'c') 

with 

[GSFC01* = r ° vv . ,. * ■ < 21 > 



where q is the on-shell momentum satisfying the equation: 



^M 2(LJ+1) A<? J BW) 2 + (A J BW ) 

9bw/ V ? 2 +(A 3 BW ) 2 



W = v/<+? + V m Lw + ? 2 > (22) 

and ?tir bw = 770 MeV; q = q^ w when = Mg W . The integer L J is the lowest allowed 
orbital angular momentum between the pion and a vector boson for a given M* at rest. It 
is noted that the BW M* propagator in Eq. (121 p is diagonal with respect to the index j, 
while the M* propagator for the UM given in Eq. ( 1121) can have off-diagonal component, 
connecting different M*s belonging to the same quantum number. Furthermore, whereas the 
BW M* propagator in Eq. (12 ip is purely phenomenological, the M* propagator for the UM 
needs to be given by Eq. ffl~2]) as a consequence of three-body unitarity. The IM defined by 
Eqs. (T2Q|) and ( l2Tj) will be used to extract the resonance parameters from fitting the "data" 
generated from our UM defined in Sec. Ill BL 
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III. RESULTS 



Our objective in this paper is to investigate the importance of the three-particle unitarity 
in determining the excited meson M* properties from fitting the Dalitz-plot distribution 
data of 3n states for the_7p — > M*n — > 7i + 7i + 7i~n reaction. We do not make an attempt to 
analyze the CLAS data |4J for this reaction here. Thus it is sufficient to consider the VDM 
production mechanism illustrated in Fig. EJ^b). Then we can set up a model by fixing the 
parameters associated with M* states using partial widths predicted by the 3 Po model of 
Ref. @. This will be explained in Sec. HTTA1 

Once the M* parameters are fixed, we can perform the calculations using the formula 
presented in previous sections since all parameters needed to calculate the irR — > ttR' ampli- 
tudes [Eq. ([8])] and nR propagator [Eq. ffTO]) ] have been determined in our previous work ((jj]. 
The Dalitz plots for the 7r + 7r + 7r~ are calculated using Eq. fl2]), and are presented in Sec. IIIIBI 
In Sec. IIII C\ we describe how the generated Dalitz plots are used as the data to determine 
the parameters of the IM described in Sec. Ill CI 

In Sec. IIII D\ we examine the differences between the resonance parameters extracted 
with our UM and those from the fitted IM. Their differences will indicate the importance 
of the three-particle unitarity in determining the excited meson M* properties from fitting 
the Dalitz-plot distribution data of the 3n states for the 7p —> 7r + 7r + 7r~n reaction. 



A. Determination of the M* parameters 

In our UM, we consider the irR partial waves that are found to be necessary to fit 
the CLAS data for 1.0 < W < 1.36 GeV fl. We thus have one or two bare M* states 
for four partial waves: J PC = 1++ M1230), Oi(1700)], 2++ [a 2 (1320), a 2 (1700)], 2~+ 
[7r 2 (1670), 7r 2 (1800)], 1~ + [7T!(1600)]. We assume that our bare M* states can be identified 
with excited meson states of the qq excitation type listed in Ref. jsj, and that their bare 
M* — > ttR couplings are fixed so that the partial decay widths predicted by the 3 Pq model jsj 
are reproduced; we use the formula given in Appendix I of Ref. [ll[ to calculate the partial 
widths. We also assume that the daughter ttR states have the lowest allowed orbital angular 
momenta. The bare masses M^* are also identified with the excited meson masses listed 
in Ref. jsj. The only exception is the 717(1600) that is speculated to be a hybrid state. In 
this investigation we use the mass and partial widths for 77(1600) given in Ref. [To|. For 
simplicity, we set all M* — > nR cutoffs to A( 7rR ) ; A /» = 1 GeV. With the above specifications, 
the parameters for our UM are fixed. 



B. Calculations of Dalitz plots 

We use Eq. (T5]) to generate the Dalitz-plot distributions of 7r + 7r + 7r~ for the 7p — > 7r + 7r + 7r~n 
reaction at the photon energy E 1 = 5 GeV and the momentum-transfer t = —0.4 GeV 2 . 
This is the kinematics considered in the CLAS analysis Q. 

We next need to specify the Euler angles [fouler in Eq. (j^J)] that define the orientation of 
the three-pion plane in their center of mass system. Experimentally, it would be preferred 
to choose an orientation where the three pions have less chance to interact with the final 
nucleon. Thus we choose the Euler angle /3 such that the three-pion plane is perpendicular 
to the direction of the final nucleon. Because the cross section does not depend on a, we 
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FIG. 5. (color online) ^-dependence of integrated Dalitz-plot distribution of the three pions for 
the 7p — > M*n — > 7r + 7r + 7r~n reaction, obtained from Eq. ([2]) by integrating over the kinematical 
variables m 2 ab , m| c and 7. 

set a = 0. The remaining Euler angle 7 gives the rotation of the three pions around their 
CM on the plane specified by a and (3. We calculate Dalitz plots by varying 7 in the range 
< 7 < 2tt. 

To see the contributions of the considered resonances to the generated data, we show in 
Fig. [5] the cross sections calculated from Eq. (j2J) by integrating over Tt^ 6 ,m^ c and 7. We see 
a broad bump at W ~ 0.95 GeV due to ai(1230), a highest peak at W ~ 1.25 GeV due to 
a 2 (1320). The second highest peak at W ~ 1.65 GeV is due to ai(1700) and a 2 (1700). The 
gap at W ~ 1.7 GeV is due to an interference between 7T2(1670) and 7r 2 (1800). The exotic 
7Ti(1600) is clearly visible as a small spike at W ~ 1.6 GeV. Clearly it is a highly non-trivial 
task to fit the generated Dalitz-plot distribution data, in particular, in the region 1.6 GeV 
< W < 1.8 GeV where the contributions from several resonances overlap strongly. 



With the kinematics described above, the Dalitz-plot distribution of the three pions are 
generated with Eq. (f2]) by varying the variables W, 7 (one of the Euler angles), m 2 ab and m 2 c . 
We calculate the cross sections at kinematical points that are uniformly distributed in the 
space spanned by these kinematical variables. The calculated cross sections at these points 
are regarded as the 'data' in determining the parameters of the IM defined in Sec. Ill C[ 

As examples, we show in the left sides of Figs. [6] and [7| the generated Dalitz-plot distri- 
butions in the m^_ 7r+ -m^. +7r+ plane at W — 1 GeV and W = 1.76 GeV, respectively. The 
Dalitz-plot distributions largely depend on 7 both in shape and magnitude. Thus in these 
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FIG. 6. (color online) (left) Dalitz-plot distribution of 3ir [Eq. (j2J); unit: /xb/GeV 7 sr rad] from 
our UM for the 7£> — > 7r + 7r + 7r _ ra reaction at IF = 1 GeV near ai(1230) peak, cos /3 = —0.37, and 
7 = 90°. (right) One slice cut of the Dalitz plot on the left at m 2 „ + = 0.51 GeV 2 . 
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FIG. 7. (color online) Dalitz-plot distribution of 3ir at W = 1.76 GeV near 7^(1800) peaks, 
cos/3 = —0.77, and 7 = 173°. On the right, one slice cut of the Dalitz plot on the left is at 



m 



1.23 GeV 2 . The other features are the same as those in Fig. [6l 



figures, we choose 7 where the Dalitz plot has the highest peak at the value of W. 

In the right-hand sides of Figs. Eland the solid curves (unitary model) are the distribu- 
tion at a fixed m 2 + _ of the x-axis of the Dalitz plot on the left. The dotted curves (without 
Z) are obtained by turning off the Z-diagram mechanism in our unitary calculations. Clearly, 
the effect of the Z-diagram mechanism, which is the necessary consequence of three-particle 
unitarity condition, is important, especially in Fig. |6] where ai(1230) dominates. 

For fitting the data with ^-minimization, we need to assign an error to each data point. 
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TABLE I. Bare masses Mj^* (MeV) as well as couplings C^. r li\ m * (dimensionless) and cutoffs 
Kr,m* (MeV) of Eq. CED, for the M* = vr 2 (1670), vr 2 (1800), and 7r'i(1600) from the UM and those 
from the IM are compared. Rf 1 denotes the i-th bare R state with the spin L and the isospin 
/. I denotes the orbital angular momentum between Rf 1 and it. See Table I of Ref. 0] for the 



properties of Rf 1 
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1000. 


1719. 


1000 1116. 


C(- K Rl 1 ) 1 ,M* 




- 2.29 - 6.66i 




2.64 + 0.30i 


0.07 + 0.16? 


C (ttR\ 1 ) s ,M* 




- 0.02 - 0.02i 




- 0.01 + O.OOi 




^(wRl, 1 ),M* 




884. 




1007. 


525. 


C (nRf>)o,M* 


8.39 


15.09 + 1.29i 


9.29 


- 3.66 - 0.55i 




C (nRf) 2 ,M* 




0.11 + 0.02i 




0.00 + O.OOi 




A (7rRf ),M* 


1000. 


817. 


1000. 


1077. 





We use the same error for all data points belonging to the same W . For a given W, we use 
5% of the highest value in the Dalitz plot distributions as the error. 

We adjust the BW parameters Mbw, Tbw and Abw m Eq. (121~|) . and all M* — > nR 
couplings constants C^r) l ,m* and cutoffs ^(■ k r) l ,m* t° fit the Dalitz plots generated with 
our UM. First we tried fitting with real M* — > irR couplings, but no satisfactory fit is 
obtained. Thus we allow M* — > ttR couplings to be complex, as usually done in the IM 
analysis. To get high precision fits, we find that the IM needs to include more partial waves 
in the M* — > ttR transition. This can be seen in Table [J in comparing the parameters of 
the starting UM and IM. The comparisons of the parameters for the other M* states are 
similar and therefore need not to be given here. (Note that the parameters of the UM are 
from selecting only few partial waves which have information from the 3 Pq model. So the 
Table [I] should not be used in comparing the merits of each model. Namely, if we start from 
the data generated from the IM with few partial waves, then the UM fits will also need more 
partial waves. ) 

We find that even with much more partial waves in the fit, as seen in Table HJ the IM still 
can not fit data well. The high precision fits are obtained only when we follow the previous 
IM analysis [IJ to add at each W a flat and non-interfering background to the Dalitz-plot 
distributions calculated from the IM. In Fig. [SJ we show the background contribution to the 
integrated Dalitz-plot distribution [m^ +7r _, ?n^+ 7r + and 7 are integrated over from Eq. (T2])] 
which can be compared to Fig. [5j The background contribution is highly TU-dependent. 
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FIG. 8. (color online) The contribution from the constant background cross sections to the inte- 
grated Dalitz plot obtained from Eq. ([2]) by integrating over m^. +7r _, m^. +7r+ and 7. 



TABLE II. Pole positions M% M from the UM and M 1 ^ 1 from the IM are compared. ' are 

the BW masses and widths (Mbw, — Fbw/2) from the IM. All are in unit of MeV. 



M* 


1V1 M* 


M™ 


M™ 


M IM(BW) 
R 


ai 


1230 


(913, -69) 


(940,-64) 


(1111,-600) 






(1443, -342) 


(1201,-212) 


(1391,-389) 




1700 


(1658,-53) 


(1672,-59) 


(1676,-59) 


r/ 2 


1320 


(1263, -21) 


(1262,-22) 


(1267,-24) 




1700 


(1652,-38) 


(1657,-48) 


(1668,-52) 


7T2 


1670 


(1786,-228) 


(1701,-220) 


(1815,-283) 




1800 


(1722,-26) 


(1724,-34) 


(1726,-34) 


7Tl 


1600 


(1599, -4) 


(1599,-4) 


(1599,-4) 



The largest contribution relative to the full Dalitz-plot distribution (the background plus 
the IM) is at W = 1140 MeV by 37%. Below the a 2 (1320) [a 2 (1700)] peak, the contribution 
is 22% [29%]. In some regions, on the other hand, the background contribution is almost 
zero as can be seen in the figure. 

A good quality of fit has been obtained with the IM. The IM gives Dalitz-plot distributions 
that are not distinguishable from the left-hand-side of Figs. |6] and [7] of the UM. The good 
fits can be seen from comparing the solid (unitary model) and dashed (isobar fit) curves in 
the right-hand-side of the same figures. 

D. M* Parameters 

We apply the analytic continuation method developed in Ref. [l5[ to extract the resonance 
pole positions from searching for the solutions of det[G~jJ*(E)] = of Eq. (TT21 for our UM 
and [G^Y(E)]^ 1 = of Eq. fl2D for the IM. The extracted poles M R for the UM and M^ M 
for IM are compared in Table ILT1 For the UM, we also list their bare masses M^*. For a\, 
we see that two bare states evolve into three resonance poles. The situation is similar to P n 
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nucleon resonances reported in Ref. 16|. In the last column, we also list the BW positions 
(M BW ,-T BW /2). 

The results shown in Table [Til are similar to what we have observed in other analyses: if 
the data are fitted equally well, the extracted resonance positions are rather insensitive to 
the parametrization of the amplitudes as far as the singularities of the scattering amplitudes 
like branch points are far from the resonance pole. Here we find a significant difference 
between UM and IM in two poles associated with second a\ (second row in Table HTj) . The 
pole position of a\ at 1443 — z342 MeV of the UM is close to the branch points at 1387 — 101z 
MeV for the tt-/ 2 (1270) channel and 1487- 167z MeV for the tt-/ (1370) channel. Thus the 
resonance shape of the amplitude at real W are distorted by those singularity at the complex 
energies. On the other hand, BW model has only a single complex energy surface, and has 
no branch cut. We have observed [l6[ the similar situation in the study of Roper nucleon 
resonance TV* (1440). 

For a2, we see from Fig. [5] that two resonances for this partial waves are well separated. 
Furthermore, a2(1320) is a very pronounced and isolated resonance, like the A(1232) reso- 
nance. Thus it is not surprising to see that the resonance positions for a2(1320) from two 
models are almost identical and the BW value is also very close to the pole position. 

For 7T2, we assign the wide resonance as 7r 2 (1670) and the narrow one as 7r 2 (1800) in 
Table [TTJ In UM, the mass of the wide resonance (1786 — 228i MeV) is higher than that 
of the narrow one (1722 — 26* MeV), while the order of two resonances is reversed in IM. 
Though IM reproduces well the data of Dalitz plot, the resonance positions for overlapping 
resonances are sensitive to the reaction dynamics. 

We next compare the residues of the partial-wave amplitudes of 77T — > Rir transition 
evaluated at the resonance positions. To proceed, we first note that the total amplitude of 
7V — > ttttttN' defined by Eqs. (JI])-(J7j) can be cast into the following form: 

cyclic 

c' R',c"R" T c ii ^j" TV' ,7 V 

(W) , | 7 JV > (23) 

(a'b'c 1 ) R' c"R" 



with 



T° c , R , N ,, lN {W) = E E T ^^L i G M4W)} 3ak T MLN , nN (W) . (24) 



J 



By using Eq. ffTUj) for T M * N r aN (W) and performing the partial-wave expansion, the following 
partial wave amplitude ^ r li^ l (W) can be separated from T® R , N , lN {W): 



-T M * ka ,( P >ir) L .(p P >) , (25) 



where Lj (Lf) is the orbital angular momentum of 771 (i?f i tt). The symbol p R Li (p p i) denotes 
the on-shell momentum of Rf 1 (p') that satisfies W = M R = E r li{v r li) + E n (p R Li) where 

i i i 

the bare mass of R\ l in our model is used to calculate the energy E r li{j> r li). 
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The residue at each pole position M R is then defined by a contour integration along a 
closed path (Cm r ) around the pole position Mr. We further multiply it by phase-space 
factors to define 



B a ' MR (M Li (RFir) Lf ) = -Jp,AM R ) Pr^AMr) <f dW T^ )l ^ [r ^ )l (W) 06) 

Jc Mr 1 

where the phase space factors are 



/ w n PrEr{Pr)E*{Pr) 
PrA w > = 71 , 

PyAW) = • ( 27 ) 

For the IM, we can cast Eq. (120]) into the same form of Eqs. (123|) -( |25|) . It is straightforward 
to see that the resulting 77T — > ttR amplitude is 



T^ L ]^ R L l7v)Lf (W) =Y. r (^h r M L (PRfi) 



y \<%Zm\u. [j;T, ] l^^ L , v , K Mj . (28) 

where Gf}¥(W) is given in Eq. ( 12T|) . Its residue can be measured by using Eq. ( 126]) with 
the replacement T^^W) T^^W). 

Following the usual procedure, the residue of the BW parametrization at W = M BW is 
then obtained from Eq. (1281) by replacing G^(W) with 1/2. Multiplying the same phase 
factors, we then define for the BW: 



where the on-shell momenta are for the BW mass (Mbw) taken as the total energy. 

We have determined the residues for M*s of UM and IM using Eq. ( |26i) . and also the 
conventional BW residues for IM using Eq. ( 129|) . In Table II I II we present only results 
which are useful in revealing the essential differences between these three residues. Because 
the overall phase is arbitrary, we choose the overall phase for IM such that the phases of 
.8(77? — > -Rpvr) for the most prominent a2(1320) are the same for all three cases in this 
table. Guided by the results of the H^-dependence of cross sections shown in Fig. [5] and the 
pole positions shown in Table [Til we focus on three different situations: (a) the resonances 
are broad such as the second ai(1230); (b) resonances in the same partial wave are narrow 
and isolated such as a2(1320) and 02 (1700); (c) resonances in the same partial wave overlap 
such as 7r 2 (1670) and 7r 2 (1800). 

We see that for the broad resonance 2nd-ai(1230), the extracted residues are drastically 
different between UM and IM. This is of course partly due to the difficulties in getting the 
same pole position for the reason discussed in this subsection. For the 02 (1320) that gives 
the very pronounced peak in Fig. [51 both the resonance positions and residues agree almost 
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TABLE III. The pole positions Mr and the weighted-residues B^jtt — > ttR^ 1 ) defined in Eq. (|26p 
from UM and IM fit are compared. Rf 1 denotes the z-th bare R state with the spin L and the 
isospin /. For Mr, they are listed as (Re(M/j), Im(Mjj)) [MeV]. The column under IM(BW) are 
from IM but using the usual BW procedure with their masses listed as Mr = (Mrw, ~ 
For the weighted-residues, we list as (|-B| [MeV], 4> [deg]) of the expression B{^ir — > -kR) = \B\e 1 ^ . 
The overall phase for IM is chosen such that the phases of B^ir — > R^tt) for the most prominent 
a2(1320) are the same as UM. 



M* 










UM 


IM (pole) 


IM (BW) 


2nd-ai(1230) 


Mr 








(1443, -342) 


(1201, -212) 


(1391,-389) 






B(^it - 


■» ttR] 1 ) 


(0,0) 


64.1, -67. 


28.4, -171. 


33.6, -150. 






B^TT - 


> ttRI 1 ) 


(0,2) 


8.1, 34. 


0.0, 180. 


0.0, 180. 


a 2 (1320) 


Mr 








(1263,-21) 


(1262,-22) 


(1267, -24) 






B(j7r - 


+ irRl 1 ) 


(2,2) 


6.2, 171. 


6.2, 171. 


6.7, 171. 


a 2 (1700) 


Mr 








(1652,-38) 


(1657, -48) 


(1668, -52) 






B^TT - 


+ nR\ l ) 


(2,2) 


9.4, 147. 


6.5, 139. 


7.1, 140. 


vr 2 (1670) 


Mr 








(1785,-229) 


(1701,-220) 


(1815,-283) 






B(-~f7T - 


+ irR™) 


(1,2) 


2.0, -123. 


4.5, 129. 


5.4, 176. 






B^TT - 


+ ttRI 1 ) 


(1,1) 


31.8, -134. 


37.6, 167. 


38.6, 168. 






B(jTT - 




(1,3) 


1.9, -101. 


0.6, 153. 


0.7, 176. 


7T 2 (1800) 


Mr 








(1722,-26) 


(1724, -34) 


(1726, -34) 






B^TT - 


+ TTii? ) 


(1,2) 


0.6, 73. 


0.6, -26. 


0.6, -35. 






B^ir - 




(1,1) 


4.5, 4. 


4.3, -7. 


4.3, -21. 






B^ir - 




(1,3) 


0.2, 14. 


0.0, 180. 


0.0, 180. 



perfectly between UM and IM. This is not surprising since it is similar to the situation 
of the well known A(1232) resonance in the irN scattering. For the a 2 (1700), there is 
some difference in resonance positions and residues. Since the residues from UM and IM 
are from loop integration around the pole position as given in Eq. fl26|) . their differences 
originate from the differences in their amplitudes near the pole position. This is illustrated 
in Fig.[9]for this resonance. We see that the amplitudes from UM and IM are rather different 
near their pole positions and hence lead to the differences in residues which are calculated 
from loop integration of the amplitude around the pole position. This indicates that the 
resonance properties extracted from data are not independent of analysis method, and thus 
it is essential to have a parametrization of the amplitude with theoretical constraint. Finally, 
for 7r 2 (1670) and 7r 2 (1800), it is not surprising to see that the phases of their residues from 
UM and IM are very different since two resonances are overlapping. 

IV. SUMMARY 

We applied the unitary coupled-channels model developed in Ref. j^j to investigate the 
issues concerning the extraction of meson resonances from the three-pions photoproduction 
reaction on the nucleon. Our aim here is to examine the importance of the three-body 
unitarity, which is not accounted for rigorously in the commonly used IM analyses for the 
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FIG. 9. (color online) Three-dimensional and contour representations of the real part of 
(W-M R )T? . , Rll . (W) [10~ 6 MeV" 1 / 2 ] for a 2 (1700) from the UM (left), and from the isobar 
plus BW model (right). The resonant amplitude, T?y^, Rll ^(W), is defined in Eq. (l25j) . The 
boundary of the W-region shown in the figure is given by the circle with the center at Mr and 
~10 MeV radius. 

resonance extraction. This has been done by comparing the resonance parameters extracted 
with our UM and an IM both of which reproduce the same Dalitz plots over relevant kine- 
matical region. We also compare the resonance parameters with the usual BW parameters 
of the same IM. 

We found that the good IM fits to the Dalitz plots generated from the UM can be achieved 
only when the M* — > irR coupling are allowed to become complex and a flat background is 
added at each W. The resonance positions from the two models agree well, except for the 
resonance poles whose tails are partly blocked by branch cuts on the Riemann surface from 
reaching at the physical real W axis and for the overlapping resonances. The residues of the 
resonant amplitudes extracted from the two models and those from the usual BW procedure 
agree well only for the isolated resonances with narrow widths. Most of the extracted residues 
for overlapping resonances could be drastically different. Our results suggest that even with 
high precision data, the resonance extraction should be based on models within which the 
amplitude parametrization is constrained by three-particle unitarity condition. 
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